Structure of Amplitude Correlations in open chaotic Systems 
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The Verbaarschot-Weidenmiiller-Zirnbauer model is believed to correctly represent the correla- 
tions of two S-matrix elements for an open quantum chaotic system, but the solution has considerable 
complexity and is presently only accessed numerically. Here a procedure is developed to deduce its 
features over the full range of the parameter space in a transparent and simple analytical form 
preserving accuracy to a considerable degree. The bulk of the VWZ correlations are described by the 
Gorin-Seligman expression for the 2-amplitude correlations of the Ericson-Gorin-Seligman model. 
The structure of the remaining correction factors for correlation functions is discussed with special 
emphasis of the role of the level correlation hole both for inelastic and elastic correlations. 

PACS numbers: 05.45.Mt,24.60.Dr,24.60.-k 



I. INTRODUCTION 

The Verbaarschot-Weidenmiiller-Zirnbauer (VWZ) 
model for chaotic scattering for an open system |l| 
is generally believed to be a realistic example of such 
processes 0, It is based on a Random Matrix Model 
for which the random Hamiltonian is dynamically cou- 
pled to the external channels. For the correlations of two 
S-matrix elements it gives an explicit general solution, 
which unfortunately is very complex. In spite of the im- 
portance of this result, the dominant physical features 
are mostly not transparent. Hence they cannot be reli- 
ably exploited in simplified form for application to higher 
order correlation functions such as cross-section corre- 
lations for which the corresponding VWZ solutions are 
not known presently. The VWZ amplitude correlations 
can presently be evaluated only using numerical methods 
with specified values for the numerous parameters [U, H| , 
but for special cases |6j. Approximate expressions exist 
for small correlation times, but the convergence radius of 
the expansion is small @, HJ • 

The properties of higher order correlation functions in 
the VWZ model are largely unaccessible with present 
approaches and their properties are unexplored with few 
exceptions. The cross-section correlations (4-point am- 
plitude correlations) have been obtained for a vanishing 
correlation time [9l4ll|. Their behavior is however only 
accessed numerically. A realistic study of these corre- 
lation functions requires then a simplified approach in 
which the properties of the amplitude correlations closely 
mirror those of the VWZ model. In the following we 
will demonstrate that the model due to Ericson-Gorin- 
Seligman (EGS) d, [HJ gives a good general represen- 
tation of the VWZ solution for amplitude correlations, 
but with some characteristic deviations. Its explicit so- 
lution is in some respects more general than the one of 
the VWZ model and has been used in a study to be 
published to explore the structure of cross-section corre- 
lations both for near-chaotic and for chaotic scattering 
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from open systems [13j. 

The EGS model is inspired by reaction theory such as 
in Ref. [Tij . It approximates the S-matrix by a sum of 
complex poles, the positions of which are assumed not to 
be influenced by the coupling to the external channels. 
In its main version these poles have partial width ampli- 
tudes with a normal (Gaussian) distribution; the pole po- 
sitions are assumed to have the same spacing distribution 
as for the Random Matrix Model [15j . i. e., generated by 
matrix elements with a normal distribution. Although 
the EGS model implements a number of features of the 
VWZ one, it differs in particular by not accounting fully 
for the dynamical feed-back of the open channels on the 
spacing distribution 16]. In addition, the price for the 
solvability of the EGS model is that it violates unitar- 
ity sum-rules, which potentially may produce unrealistic 
properties. It is shown however in Ref. [l3[ that these 
restrictions are of unimportant in most situations; this 
conclusion is concretely vindicated by the results of the 
present article. As a check on the effectiveness of this ap- 
proach the Gorin-Seligman (GS) results for the 2-point 
correlations of the EGS modelQ have been numerically 
compared to the ones in the VWZ modelQ. The find- 
ing is that this model with normally distributed partial 
width amplitudes describes the gross features of the VWZ 
model surprisingly well for a variety of conditions. This 
suggests that the structure of the GS result might be 
used as a guide to develop an accurate and transparent 
analytical approximation to the VWZ results. This is the 
goal of the present article. It gives a valuable insight in 
addition to the accurate numerical results for the 2-point 
function. This permits one to identify explicitly the ori- 
gin and nature of the differences of the VWZ result and 
the GS one. The understanding of this point gives con- 
fidence that the EGS model gives a reliable approximate 
description of the chaotic cross-section correlations [HI . 
The strategy, therefore, will be to extract the main varia- 
tion of the key VWZ result with the GS solution as guide 
to identify appropriate expansion parameters for the de- 
viations. Using this information a natural procedure is 
developed for obtaining explicit closed expressions for the 
additional factors over the entire parameter space. 
In Sec. |H]the standard VWZ 2-point function (8.10) in 
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Ref. [1] is converted to the time representation and ex- 
pressed more conveniently in terms of the time variable r. 
This representation is technically easier to handle than 
the energy one, since it gives products and not folding 
integrals. In addition, it is immediately apparent that 
there are no singularities in the integrals. The basic cor- 
relation functions and their notation are defined here. 

In Sec. IIIII some basic identities associated with VWZ 
integral are given as well as their relations to the level 
spacing correlations. This section also defines notations 
for the weighted functions generated by the integrations. 



Sec. IIVI describes the construction of the small expan- 
sion parameters, the basic expansion and the expansion 
of the factors coupling to the observed channels referred 
to as the 'leads'. It gives in addition the structure of the 
general expansion term. 

Sec. [V] gives the explicit expressions for small and large 
r as well as a convenient interpolation connecting the two 
regions. 

Sec. |VT] is devoted to a discussion of the characteris- 
tic features of the amplitude correlation functions in the 
light of the results. It also discusses possible generaliza- 
tions and some practical consequences. 



II. THE VWZ CORRELATION FUNCTION 



The key result in Ref. [1] is Eq. (8.10), which gives the correlation function for two conjugate S-matrix amplitudes 
S(E(1)) and S*(E(2)) at energies E(l, 2) differing by e = E{\)—E(2). The system is assumed to have no secular energy 
dependence with [E(l) + E(2)}/2 and threshold effects in the different channels are neglected. The observed channels 
(the 'leads') are denoted by a, 6, c, d while open channels without restrictions are denoted by e with corresponding 
transmission coefficients T e . Ensemble averages are denoted by < >. The energy scale has units such that the 
average level spacing d = 1. 

The general VWZ correlation function in the energy representation is given by the triple integral 



C ab - M {e) =< S ab (E(l))S* cd (E(2)) > - < S ab >< S* d >-- 



(1) 



dXo 



dX- 



(l-A)A|Ai-A 2 | 



o [(1 + A 1 )A 1 (1 + A 2 )A 2 ] 1/2 (A + A 1 )2(A- 
exp [-(wre)(Ai + A 2 + 2 A)] ■ 0. - L X) 



{ 



SabScd < S aa >< S* > T a T ( 



(l + T e A 1 )i/2(i +Te A 2 )i/2 

Ai A 2 2A 



{SacSbd + SadSbc) T a Tb 



l + T a Ar 
Ai(l + Ai) 



1 + T a A 2 1 - T Q A 
A 2 (l + A 2 ) 



Ai 



A, 



2A 



T c Ai 1 + T C A 2 
2A(1 - A) 



1-T C X 



(l+TaAOa + TjAi) (l + T a A 2 )(l+T 6 A 2 ) (1 - T a A)(l - T b X) 



The normalization is determined by the triple integral of a weight factor and gives no singular contribution from the 
region near the origin 



dXi 



dX 2 



dX 



(l-A)A|Ax - A 2 | 



o [(1 + Ai)Ai(l + A 2 )A 2 ] 1/2 (A + Ar) 2 (A + A 2 ) 2 



= 1 



(2) 



While these expressions present nicely in these variables, it is not the most convenient representation for applications. 
The equivalent correlation function in the time representation is obtained by taking the Fourier transform using that 
/_ de exp[27rixe] = 5(x), where the energy is measured in units of the average level spacing d. The time r > is 
measured in standard units 2ir/d as , e. g., in Ref. [12], Eq. (8). This means in particular that the condition for 
overlapping levels T > d corresponds to ^ e T e > 2n . It is convenient to make the variable substitutions Ai >2 = 2txi )2 
and X = tx . These transformations give immediately the Fourier transform as 



C ab -cd{r) 
dxi 

5 [1 - (xi + x 2 



de exp(27rier)C £lb;cd (e) 



(3) 



dxo 



min(l,l/r) 



dx 



(1 - TX) 



x\x\ — x 2 \ 



"in a 



(1 + 2tx 1 ) 1 / 2 (1 + 2rx 2 ) 1 /2 (xix 2 y/ 2 {x + 2 Xl ) 2 {x + 2x 2 ) 2 
(1 - rT e x) 



Sabred < S aa >< S* > 2rT a T ( 



(1 + 2ra; 1 T e ) 1 /2(i + 2ra; 2 T e ) 1 /2 
Xi x 2 



1 + 2rT a a;i 1 + 2rT a x 2 1 — rT a x J \ 1 + 2rT c a;i 1 + 2tT c x 2 1 — tT c x 



Xi 



Xn 



+ 
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(SacSbd + Sadfibc) T a T b 



xi(1 + 2txi) x 2 (1 + 2tx 2 ) x(l — tx) 



(l + 2 T T a x 1 )(l + 2TT b x 1 ) {l + 2TT a x 2 ){l + 2TT b x 2 ) (1 - tT s)(1 - rT b x) 



The expression ([3]) is the basis for the following. It is convenient to work with the correlation functions C\ a . cc (T) and 
Cab a bi T ) defined from this expression: 

C a b-cd( T ) = S abS c dCa a . cc (T) + (S ac S b d + 5 ad 5 b c) Cl b . ab {r) (4) 

There are 3 basic situations: elastic autocorrelations , correlations of two different elastic amplitudes and inelastic 
autocorrelations. Other correlations between two S-matrix amplitudes vanish for the VWZ 2-point function. The 
relevant cases are explicitly 

Caa:aa( T ) = CL-,aa( T ) + 2 CL-,aa( T ) elastic autocorrelations (5) 

C a a:cc{ T ) = ^l.a-cc{ T ) elastic correlations with a/c (6) 

C a b;ab(T~) — C a b;ba{T~) — ^ab-ab( T ) inelastic autocorrelations with a ^ b (7) 



III. BASIC IDENTITIES 



The following identities are crucial in the expansions for small and moderate r. They result in the limit of vanishing 
transmission coefficients and are central in deriving non-trivial approximations in the following. Exact integrals are 
obtained from Eq. ([3]), where the index r on a quantity means that it is an averaged integral weighted by the parts 
of the integrands that are independent of the transmission coefficients as in Eq. ([5| below. In the limit of vanishing 
transmission coefficients the correlation function C\ (r) defines the following function 



J(r) 



(8) 



dx\ 



dxo 



mm(l,l/r) 



dx 



(1 - TX) 



x\xi - x 2 \ 



(1 + 2tx 1 ) 1 / 2 (1 + 2tx 2 ) 1 ' 2 (xixa) 1 / 2 ^ + 2 Xl ) 2 (x + 2x 2 f 



x 5[1 - (xx + x 2 + x)] 



Here I(r — 0) = 1. One notes that the time-dependent factor formally appears as an additional channel e with T e = 1. 
We will in the following denote by [f(x±, x 2 , x)] t I(t) the integral with the same time-dependent integrand weighted 
by f(xi, x 2 ,x). In Eq. (jHJ) this definition of I(t) corresponds to the integration of the integrand multiplied by 1 with 
[1] T = 1. " 



2tI(t) = 1 -6i(t) 



(9) 



The last term in this expression for J(r) is exactly the Dyson level correlation function bi (r) [Tj] , while the first term 
1 in Eq. ([9]) comes from the self-correlation of a narrow level with itself. Here b\{r) = J dsY^i^s) exp(27r«Ts) of the 
corresponding level correlation in the energy representation. It has the explicit form [15| 



For < r < 1 6i(r) = 1 - 2r + rln[l + 2r] ~ 1 - 2r+ 2r 2 - 2r 3 + 8r 4 /3 - 4t 5 



For 



t > 1 b x (r) = -1 + rln 



2r + 1 
2t- 1 



+ (l/12)r" 



(10) 



Here &i(l) — 0.0986.. and it vanishes rapidly for larger r. 
This level correlation function is an explicit result in the 
VWZ model found by Efetov and Verbaarschot [H [l|[ . 

A similar, but simpler, relation is obtained similarly 
from the correlation function C 2 (t) in the limit of all 
T e = 0. In this case the origin is also a level correlation, 
but here only the self-correlation of the levels contributes 



and gives only the first term in Eq. ([9]). 

[1 + 2t{x\ +x 2 2 - (l/2)x 2 )] T /(r) = 1 (11) 

It is convenient to regroup the terms in Eq. ([lip using 
x\ + x 2 + x = 1 so as to display the parameters x\X 2 and 
x which will serve as small expansion parameters in the 
following: 



x\ + x\ - {\j2)x 2 = 1-2 (x 1 x 2 +x - (l/4)x 2 ) 



(12) 
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With the expression (j9j) for 7(r) in terms of the Dyson function 

At [xix 2 +x-{l/A)x 2 ] T I(t) = (1 + 2r)/(r) - 1 



l-(l + 2r)6 1 (r) 



2r 

(13) 



r(l-r + 2r 2 /3 + 0(r 3 )) (for r < 1) ; ^ - ^ (/or r > 1) 

This combination is used later in the perturbative expansion terms. 

IV. STRATEGY FOR EXPANSIONS 

Our method is illustrated by the inelastic case assuming at first for simplicity that the transmission coefficients in 
the leads T a ^^ c are small. In the case of the GS result, the dependence on the transmission coefficients in Eq. (IA1|) 
is then governed by the factor n e (l + 2rT e )- 1 / 2 . A related factor J] e (l - txT c )(1 + 2tx 1 T e )- 1 / 2 (\ + 2rx 2 T e )~ 1 / 2 
appears in the VWZ model as the product contribution to the integrand of Eq. which approximately gives such 
a factor for small and moderate T e . This suggests that one should avoid power expansions based on 2rT e < 1 as a 
small parameter, since such expansions mostly have a narrow radius of convergence for the GS case. A far better 
expansion is obtained by the following procedure. 

Since x\ + x 2 = 1 — x and with the notation T' e = T e /(1 + 2rT e ) one can write after breaking out the factor 
(l + 2rT e )- 1 /2 : 

(l- T T e x) (l-rT e x) iM , o_^-i/2 



(1 + 2rxiT e ) 1 /2(i + 2rx 2 T e ) 1 / 2 (1 - 2rT' e x + Ax 1 x 2 T 2 T e T' e ) x l 2 



x (l + 2TT e y l,z (14) 



The first term of this product is then exponentiated. The individual contributions in the exponent are expanded 
assuming sufficiently small expansion parameters 

(1 - rT e x) = 

(l + 2Tx 1 T e ) 1 / 2 (l + 2Tx 2 T e ) 1 / 2 1 ' 

exp |^ [ln(l - rxT e ) - (1/2) In (l - 2rT' e x + A Xl x 2 T 2 T e T' e )] j x + 2rT e ) _1//2 ~ 
exp{-2r 2 T e T' e [x x x 2 + x - (l/4)a; 2 )} x JJ(1 + 2rT e )- 1/2 

e e 

In this expansion the linear terms in x\x 2 and x regroup naturally with the term in x 2 for small r. For r > 1 this is 
no longer the case, but terms of second order and higher in these expansion quantities are negligible in the relevant 
regions of integration as compared to the linear terms. 

Note that the expansion is unrelated to the issue of whether the resonances are overlapping. 

Corresponding approximations can be made for the lead correction factors g ( j i (r) which according to Eq. Q are 

it \ ( Xl x 2 x \ ( xi x 2 x \ 

g 6l {{x 1 ,x 2l x;T) = + — —-= + — x — + — —-= — + — (16 

\l + 2rJ a a;i l + 2Tl a x 2 \~T± a xJ \1 + 2t1 c xi 1 + 2ti c x 2 1 — t1 c xJ 



gc2{{x 1 ,x 2 ,x\r) 



xi{\ + 2tx\) 



x 2 (l + 2tx 2 ) 



x(l 



(1 + 2tT q xi)(1 + 2tT 6 xi) (1 + 2rT Q x 2 )(l + 2 T T b x 2 ) (1 - rT a x)(l - rT b x) 



(17) 



Assuming as previously that no single decay channel dominates contribution from the lead channels and extracting 
an overall factor (1 + 2TT a j ) c )~ 1 one finds a qualitatively different contribution, which now contains terms which are 



The following expansion assumes that 2x\X2rT e < 1, typically. 
Equation II13I I suggests the characteristic value of x\X2 is of order 
0(4t)~ 1 ) or less in typical integrals for r > 1. One can verify 



that deviations lead to small and negligible contributions to the 
correction terms even in the limit of very strong transmission 
T e ~ 1. 
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linear in rT a ^ tC which is of significance for the behavior of the correlation functions for small r. This is contrary 
to the quadratic dependence on r 2 T 2 in the sum r 2 ^ e T 2 of Eq. (fT5")) . The next order terms in the transmission 
coefficients of the leads are readily included, but do not qualitatively change the results. In accordance with the 
previous assumption that no channel dominates and for simplicity only the linear contribution from the leads is 
retained here. The correction term to this order in the integrand becomes for the lead factor g^ (r). 



2Tg dl ((xi,x 2 ,x; r) ~ 2r(l + 2rT Q )" 1 (l + 2rT c )- 1 [l + 4r(T a + T c )( Xl x 2 +x- x 2 /4) + .. 



(18) 



For 5(5 2 (r) the leading terms in rT 0) i, are with r' = r/(l + 2r) 

g d2 ( Xl ,x 2 ,x; t) ~ (1 + 2r)(l + 2rT a )- 1 (l + 2rT b )- 1 [l - 4r'(l - (1 + r)(T a + T 6 )) (nu + x - a; 2 /4) + ...] (19) 

The factors given in Eqs. (|16H17p with Eqs. (fl~8| and (fT9|) . respectively, will be used in the following integrands. 

The relevant terms weighted in the VWZ integral therefore occur as exponents proportional to x\x 2 + x~x 2 /A char- 
acteristically scaled by a multiplying term dependent on the transmission coefficients. When the linear approximation 
to the lead correction is valid over the important regions of the integrals we can exponentiate them and combine the 
results with those of Eq. (|15[) . Introduced into the integrand of Eq. © gives 

CLcdr) - (Saa) (S* cc ) 2rT a T c (l + 2rT a )- 1 (l + 2rT c ) _1 TT(1 + 2rT e )- 1 / 2 [exp{- Kl ( Xl x 2 +x- (l/4)x 2 )}l 7(r) (20) 

L J r 

e 

with Ki = -4r(T a + T c ) + 2t 2 ^ T e T e ' (21) 

e 

and 

C^-aiW - (! + 2r)T a T b (l + 2rT a )- 1 (l + 2rT;,)~ 1 TT(1 + 2rT e )- 1 / 2 [exp{- K2 (x l2 ; 2 + rs - (l/4)x 2 )}l 7(r) (22) 

L J r 

e 

wtft k 2 = 4r' - 4r'(l + r)(T a + T 6 ) + 2r 2 ^ T e T^ (23) 

e 

This approximation is expressed here in a global form valid for all values of the time variable r. It will now be 
evaluated in detail in different regions. 



V. EXPLICIT RESULTS FOR THE APPROXIMATION 
A. Small and moderate r with 0(ki,2) < max[l,r] 
From Eqs. dU) and P^|) this gives expanding Eq. (j2"0"T) for the exponent smaller than unity 

C l aa , cc (r) ~ (S aa ) (S c * c > 2rT a T c (l + 2rT a )- 1 (l + 2tT c )~ 1 JJ(1 + 2rT' e )~ 1 / 2 x 

e 

I- r(T a + T c ) - -t ^ T e T e — 



2r 



2r 2 



The restrictions from the integration region are automatically takien into account. 
The corresponding result for C 2 b . ab {r) horn Eq. (f22|) is 

C 2 b;ab (r) ~ T a T h (l + 2rT a )- 1 (l + 2tT{,)~ 1 [](1 + 2rT e )~ 1 / 2 x 

l-(l + 2r)6 1 (r) 



1 + ^(1 + r)(T a + T b ) - ir 2 (l + 2r) £ T e T^ 



2r 2 



(24) 



(25) 
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For r > 1 the Dyson function in Eq. (|T0| is very small with 61 (r) ~ l/12)r 2 and gives a negligible contribution. 
Within the expansion region C\ a . cc (r > 1) becomes 



CL ;cc (r > 1) 



(26) 



(Saa) TO T„T C (1 + 2rT a )- 1 (l + 2TT,)" 1 TJ(1 + 2rT e )- 1 / 2 x ^1 + (r a +T C - -r^ T T^ (1 + ©(r- 1 ))^ 

-> (S aa ) (5 C * C ) T a T c (l + 2rT a )- 1 (l + 2tX' c )~ 1 TJ(1 + 2rT e )- 1 / 2 x ( 1 + T„ + T c - \j2 T * ) f orT ^°° and £ T e < 2 

e \ e / e 



The corresponding result for C ,2 b . ab (r) is 



C 2 b;ab (r > 1) ~ T a T b (l + 2rT a )- 1 (l + 2rT b )~ 1 ^(1 + 2rT e )- 1 / 2 x ^1 + I (t„ + T b - r £ T e T^) (1 + 0(0)^ (27) 
r a T b (l + 2rr Q )- 1 (l + 2rT b )- 1 [](l + 2rT e )- 1 / 2 X ( 1 + i(T« + T 6 ) - £$^T e j /or r -> 00 and ^T e < 2 

e \ e / e 



Note that the limit t — > 00 applies only to the case of non-overlapping levels. 

I 



B. Moderate and large r with 0(^1.2) > max(l,r) 

In this limit the previous expansion breaks down and 
the problem is dominated by the exponent appearing 
in Eqs. (|2T)|) and (f2"2")l . In its simplest form the lead 
terms are weak. Then the unobserved channels determine 
the correlations. The typical integrand has the factor 
n e (l + 2rT e )- 1 /2 [ ex p [-« 1(2 (x t x 2 + x - (l/4)z 2 )]] . It 
contains in addition the contribution to the integrand 
from the time-factor which formally corresponds to an 
additional channel with T e = 1 which can be expanded 
as for any T e with r' = r/(l + 2r). These two decreasing 
factors set the scale for contributions such that for the 
weakly coupled inelastic case these come effectively from 
the region for which 



(l/4)x 2 



O 



1 



TT , + T 2j- e T e Tl 



(28) 



Since x\ + X2 + x = 1 and x\,i enter symmetrically we 
can take x\ > x 2 and multiply by 2 . When the inte- 
gration variable in Eq. (f28|) is limited to small values, 
we have X\ — > 1, while x 2 and x are both effectively 
small with x 2 « l;x << 1. In this limit the integra- 

I 



tion can be carried out explicitly. Note however that for 
t > 1 we have x < 1/r, while the integration region 
for x 2 is limited only by x 2 < 1/2. The full integration 
region corresponds approximately to a long narrow rect- 
angle for large r. Equation (l28l) can be used to chose 
t such that the contributions come effectively only from 
a region for which both variables are much smaller than 
1/r. These conditions should be introduced into the ba- 
sic integration as given in Eq. (|SJ) but in addition the 
explicit time-dependence of the integrand should be ex- 
pressed correspondingly and exponentiated. 
Under these conditions the integral over x\ is immedi- 
ate. The integrand over x 2 ,x concentrates well within 
the integration region which can be extended to 00 for 
both x 2 ,x without restrictions. From Eqs. (|20[) and 
(f2"2")l this corresponds to the explicit evaluation of the 

terms exp — n\ , 2 (xix 2 + x — (l/4)a; 2 ) /(r) using that 

-I T 

the time-factor is (1-tx) (1 + 2rx 1 )" 1 / 2 (l + 2ra 2 )" 1/2 
~ (l + 2r)~ 1 / 2 exp[-2-n-' (x x x 2 + x- (l/4)a; 2 )] to a good 
approximation in this region. Consequently, as described 
in Appendix IBI we have in terms of the integral I (a, r) of 
Eqs. JHTJ and JB2} 



Ct;ecW - (Saa) (S*J 2rT a T c {\ + 2rT a )" 1 (l + 2tT c )~ 1 [](1 + 2rT e )-^ 2 x \ 



2 V2t 2 + (1 + 2t)(ki-1/2) 



1/2 



(29) 



where « a = -4r(T a + T c ) + 2t 2 ^ T e T' e 

e 

The corresponding expression for C 2 is closely similar : 

C 2 akab (r) ~ T a T b (l + 2r)T a T b (l + 2rT Q )- 1 (l + 2rT b )~ 1 TJ(1 + 2rT e ) _1//2 x I 

e 

where k 2 = 4r' - 4r'(l + r)(T a + T b ) + 2r 2 ^ T e T e ' 



2 V2t 2 + (1 + 2t)(k 2 -1/2) 



1/2 



(30) 
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In the region r > 1 these expressions simplify neglecting contributions of order r 1 in the correction terms. 

1/2 / 1 



C'cUccM - ( s aa) (S* c ) T a T c {l + 2rT a y 1 (l + 2tT c )- x JJ(1 + 2rT e )" 1/2 x (J^j ( i _ — - ^ - - - — j (31) 



C2**(r » 1) - T a T b (l + 2rT a )-\l + 2rT 6 )- J^l + 2rT e )~^ x (|) 1/2 ^ ± - ^ + ^ — — j ' (32) 
For asymptotically large r these expressions become 

3Ucc(t "►«>)) -> (S aa )(S* cc )T a T c (l +2rT a )-\l + 2tT c )- 1 ]1(1 + 2rT e )-^ x (^j^ (_-^__I^_^_y /2 (33) 



1/2 



1/2 



C a 2 fc;afe (r^^)^r a r b (l + 2rT a )- 1 (l + 2rT fc )- i n(l + 2rT e )- 1 /2 x _ — + ^ ) ' (34) 



2(T a +T b ) + E e T e 



1/2 



C. An Interpolation Formula 



The behavior of the correlation functions is most easily visualized in terms of a closed analytical expression valid 
for the entire parameter space, even at the cost of a slightly reduced accuracy. The previous expressions cover 
most regions of the time variable, but the small and large r expansions do not quite overlap. The interpolation 
between these two regions is easily obtained on interpolating by hand. These expressions depend on the combination 
-2(1 - r')(T a + T b )+r J2 e T e T' e for C 2 b;ab (r) and on a similar one, -2(T a + T c ) + r £ e T e T e ', in the case of <5* o;M (t). 
Except for a minor effect at small r produced by the linear dependence on the lead parameters T a ,f>,cj the correction 
term as compared to the GS expression is a smooth, monotonously decreasing positive function of r Yl T e Tl. For 
C* 2 b . ab (r), a simple interpolated expression is obtained by replacing the small r expansion in Eq. (|25[) which contains 

a term of the type 1 — y with a small parameter y by an expression of the type (1 + 2y) -1 ' 2 ~ 1 — y. 
This gives, over the entire region, the interpolated value 

C 2 ab;ab (r) * T a T b (l + 2rT Q )- 1 (l + 2rT b )~ 1 [](1 + 2rT e )- 1 /2 x 

e 

/ 1 , A (l-(l + 2r)6i(r)) I"' 72 
1 + 2( - r(l + r){T a + T b ) + -r 2 (l + 2r) ^ T e T e ') A _ L \ (35) 

In this form the same expression can be used both for isolated and for overlapping levels. This expression describes 
the global shape of the VWZ result very well over the entire range, including its absolute value even at large r for 
which the correlation functions become extremely small. It differs there from the more exact expressions by an overall 
factor 2/^/n = 1.128.. for T^rjj >> 1. For very large t the VWZ factor in Eq. (f3"5j) is dominated by the sum 
over those open channels for which rT e >> 1 since 2rJ2 e T e T' e — > J2e Te- 
la the region r > 1 nearly the same expression is valid also for the elastic correlation function C,J a . cc (r) both for 

isolated and for overlapping levels. The dependence of the correlation functions C ' on the level correlation function 
b\ (t) is nearly negligible in this region and is indicated in Eqs. (|35p and (I36[) only to emphasize the continuity of the 
expressions. For the region r < 1, however, the expression (|24[) should be used, however, so as to correctly incorporate 
the effects of the level correlations. One has 

CL; CC (r) * (S aa ) (S* cc ) T a T c (l + 2^)^ (1 + 2tT c )~ 1 Y[(l + 2TT e )- 1 ' 2 x 

e 

' \ \ -V 2 

1-(1 + 2t)6i(t) 



jl + 2( - 2r(l + r)(T a + T c ) + ir 2 (l + 2r) ^ T e T' e 



2r 2 



for r > 1 (36) 



Here the level correlation function 61 (r) is nearly negligible for r > 1, but it is included above to emphasize the 
generality of the expressions (|35|) and (|36[) for all values of r. The implications of these results will be discussed in 
the next section. 

I 

VI. DISCUSSION described by a closed analytical expression which is at a 



The previous results demonstrate that the chaotic am- 
plitude correlations of the VWZ model can be accurately 



8 



the same time both simple and transparent. 
To avoid misunderstandings, we remind the reader that 
the results have been obtained assuming for simplicity 
that no single channel dominates the sum of the partial 
widths nor the sum of their squares. This assumption is 
not basic and can be generalized. In particular, the case 
of a small number of channels can be obtained using the 
same technique even without this assumption. The price 
is a considerably more complicated discussion. 
The analytical representation reveals that the gross 
shape of the correlation function is dominated by a 
characteristic product factor which is identical to the 
one in Eq. (|A1[) for the amplitude correlations in the GS 
solution of the EGS model. This property stands out 
particularly well in the case of the inelastic correlation 
function Cab-ab( T )- For the EGS model it is a direct 
consequence of the assumption that the partial width 
amplitudes are random, i. e., have individually Gaussian 
distributions uncorrelated between channels. The result 
for the VWZ model demonstrates that in this case these 
amplitudes have become dynamically random to a high 
degree. This strong dependence on the width amplitudes 
and their distributions means that the correlation 
functions are particularly sensitive to any deviations 
from chaotic conditions in this sector and this will be 
so for the VWZ model as well. This is of considerable 
practical interest, since it opens the possibility of direct 
investigations of the sensitivity of systems approaching 
full chaoticity to amplitude distributions which are not 
normal. The consequences of possible modifications 
have been investigated in the case of the GS solution 
[l3[ . Such results are immediately relevant the VWZ 
case as well. 

Corresponding observations are also valid for the 
elastic correlation function C^ a . cc (r) which for r > 1 

rapidly converges to the same expression as for C^ b . ab (T) 
but for an overall factor. The elastic and inelastic 
correlation functions differ for r < 1 mainly due to 
the dominance of the level correlation function &i(t) in 
Eq. (fit)]) in C' 1 J Q . cc (r) in this region; a natural form for 
this function is the Dyson function which explicitly is 
a consequence of the VWZ model Q, but introduced 
phenomenologically " ad hoc! in the EGS model and its 
solutions. From the present approximate VWZ results 
one observes directly that its global properties are the 
important ones: the level repulsion on the scale of 
the level spacing and the detailed normalization of the 
"correlation hole" in the spacing distribution which 
it produces. These level spacing correlations rapidly 
become unimportant for larger r. In this region the 
correlation function explores predominantly the average 
structure of individual broadened levels. 
A general feature of the VWZ elastic correlation func- 
tion is the linear dependence on r near the origin. It 
is presently clear from Eq. ([2H]) that this is a conse- 
quence of the level repulsion by which a level creates an 
adjacent "hole" corresponding to exactly one missing 



level. Since the Fourier transform at r = of the 
correlation function C\ a - cc ij) is the energy integral 
over the entire range, the joint contribution of the 
level and its "hole" vanishes in this case. This effect 
has previously been well investigated in the context 
of the EGS modeljH; the present result reveals a 
more detailed, yet still simple, picture of the interplay 
of the level distribution and the transmission coefficients. 

It is frequently taken for granted in the literature 
that the correlation functions decrease exponentially 
for large r for strongly overlapping levels. The present 
analytical VWZ description shows clearly that this is 
only approximately correct as is also the case for the 
GS expression. The correct asymptotic dependence of 
the correlations in the VWZ model is an inverse power 
law r~( M / 2+2 ) with M is the number of open channels e 
in agreement with the behavior previously suggested by 
the EGS model j4[. With increasing r a power r^ 1 / 2 is 
added successively every time rT e becomes larger than 1 
for any channel. For experimental tests of the properties 
of such systems, it is then useful to consider not only 
the the main global exponential decay parameter ^2 e T e , 
but also J] e Tg which governs the onset of the deviation 
from the exponential law. 

The differences of the 2-point function for the exact 
VWZ solution to that of the EGS model appear mainly 
as 3 characteristic effects. 

1. The 2-point function of the VWZ model differs 
from the one of the simplified EGS model by an overall 
modulation function given in Eqs. (|35I36[) . It is a 
considerable experimental challenge to display this 
function explicitly owing to the high statistical precision 
required. Its effects are very interesting from a theoret- 
ical viewpoint, however. They can be investigated in 
detail both using exact numerical VWZ results, using 
our approximations as well as using numerical computer 
simulations of chaotic systems. This modification stands 
out particularly clearly in case of the inelastic correlation 
function C ,2 b . ab (r). It has a typical behavior apparent 
from the expressions in Eqs. (|24I25I) and (|29I30|) . 
For r > 1 the dominant term is modified by a factor 
(1 + 2rY^ e r e T e ') _1 / 2 , which decreases monotonically 
with r. In the asymptotic power law limit of very large 
r this factor becomes approximately (1 + TV TV) -1 / 2 
both for isolated and for overlapping levels^. This 
asymptotic limit is of little practical importance since 



2 A frequent practical situation is that of open channels which 
divide effectively into two groups, one with transmission coeffi- 
cients contributing both to T e and ^ T 2 at the relevant r 
and another one with small transmission coefficients contribut- 
ing only to T e , but with a negligible effective contribution to 
^2 e ■ These latter weak channels give a contribution indistin- 
guishable from the one for true absorption with an exponential 
contribution to the correlation function with a decay parameter 
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the characteristic constant factor is overshadowed by the 
power law variation t~( m / 2+2 \ 

2. The VWZ model has a typical suppression fac- 
tor < S aa X S cc > = (1 - T a ) 1 /2( 1 _ Tc )i/2 in the 
case of the correlation function C 1 contributing to 
clastic correlations as apparent already in the general 
expression (|T|) as well as in the detailed expressions 
(I24I29P and (|35I36[) . This factor has the consequence 
that the correlation function for 2 elastic amplitudes 
vanishes in the extreme limit when one or the other of 
the transmission coefficients reaches the limit of very 
strong transmission, i.e., for T a ^ c = 1. Similarly, in this 
same limit the elastic autocorrelation function is partly 
suppressed, but only by about 30%. The GS solution 
used here for comparison has the factor 1 instead of the 
factor < S aa >< S cc >, but the limit T afi — 1 is outside 
its range of its applicability as discussed in Ref. [l3j . 
The difference originates in the unitarity constraint, 
which is not imposed on the EGS model. For smaller 
T a:C the effect of this difference becomes minor. The GS 
expression can be improved by introducing this VWZ 
factor phenomenologically. One should however realize 
that the case of very strong transmission is of limited 
practical interest. This multiplicative factor in the 
elastic channels is the single most important difference 
between the VWZ and EGS models for chaotic ampli- 



tude correlations. 

3. For r < 1 the modulation function depends explic- 
itly on the level spacing distribution also in the inelastic 
case. In the region of very small r this inelastic factor is 
[l+T(T a + T b )-l/2T 2 J2 e T e +••]• At sma11 T this leads t0 
a minor initial increase of the modulation factor with a 
maximum returning to 1 near r = 2(T Q +T&)/ T 2 and 
thereafter joining the previously described monotonic 
decrease for r > 1. A closely related factor modifies 
the correlation function C£ . cc (r), but it is then difficult 
to disentangle from the level spacing correlation function. 

The present results provide an encouraging structure 
for the approximate description of higher order correla- 
tion functions for chaotic and nearly chaotic systems for 
which exact results are known only in special cases. The 
close correspondence the present results from the VWZ 
model to those of the EGS model suggests that the lat- 
ter should provide a reliable, albeit approximate, guide to 
this situation. This feature also suggests more generally 
that the simplified EGS model is a valid laboratory for 
exploring the sensitivity of various features which govern 
the approach to the fully chaotic situation of the VWZ 
model. In particular, the physics of the sensitivity of 
the correlation functions to symmetry violations can be 
reliably clarified using the simpler EGS model. 
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Appendix A: The 2-point Correlations in the EGS Model: Time Representation 



The EGS model gives the followin g re sults for the present case of normally distributed partial width amplitudes in 
each of the uncorrelated channels [il. Il3|. i.e., partial widths with Porter-Thomas distributions [l9T |: 

C[S a bS2b\(T) = T a Tf ) n e;a fc(r) for a ^ b (inelasUc autocorrelations) (Al) 

C[S aa S* c ](T) = T a T c {li e ,ac{T) - po^(T)II e;a (r/2)n e;c (r/2)^ (elastic correlations) 

for t > and for r < 0. 
Here 

n e;aA ..(r) = (l + 2rT Q )- 1 (l + 2rT /3 )- 1 ..[](l + 2rT e )- 1 / 2 (A2) 

e 



A-abs = ^ T" mal1 . In this case these weak channels do not con- 
tribute to terms VJ e T e T' e in Eqs. H31I32H nor to the asymptotic 
expression containing T e , nor do they effectively occur in the 



asymptotic power expansion. This observation in no way inval- 
idates the asymptotic power law. It only means that the region 
of exploration should be sufficiently large with rT e > 1 for all of 
the channels. 
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for a 7^ p. It is multiplied by 3 for a = j3 in the case a normal distribution. 

The exponential approximation to the time-variation holds if the variance of the total width obeys the inequality: 
(2tt) 2 t 2 < (r- < r >) 2 >= r 2 E e T 2 < 1 ■ 

Appendix B: The integral 7(a,r) in the large a limit 

Consider an integral of the type I{a 1 r), symmetrical in x\ and x 2 , for large values of the parameter a » max(l, r) 
to leading order neglecting corrections of order 1/a. 

I(a,r)= (Bl) 

rmin{l,l/T) x\x\ - X 2 \ 

dxx I dx 2 I dxS(xi + x 2 +x- I)- -rjTi — 7T, — yT( — — « ex P [~ a [ x i x 2 + x - (l/4)a; 2 ]] 

o Jo Jo (x 1 x 2 ) l -' 2 {x + 2xiY(x + 2x 2 ) 2 

For the case x\ = 1 — x 2 — x > x 2 ~ 1 it follows that the relevant region of integration over x 2 + x << 1 can 
then be freely extended to oo, since this gives negligible extra contributions. In this case the term (x + 2iei) 2 ~ 
4[1 — 2x 2 — x] — 4(xi — x 2 ) which is the numerator term such that the ratio to a good approximations is 1 which 
simplifies the discussion (an omitted term (x 2 + x/2) 2 is assumed small compared to 1). Denoting x 2 = y = u 2 and 
x = z = v 2 with p 2 = u 2 + v 2 and u = p cos </>; v = p sin 4> gives on exponentiating the contributions associated with 
x%, 

2 f°° f°° z 
I(a,r)~-/ dy dz exp [-(a - 1/2) [y + z}} = 

4 Jo Jo y 1 ' 2 {z + 2y) 2 

r-oo p7r/2 ■ 3 ± / \ 1/2 

2 l ""I ^ (l + cos^) 2 eXp[ - (fl - 1/2) ^^ KW^l)) ^ a »^(l,r) (B2) 

where the the last step is obtained using J °° dtexp(—kt 2 ) — (ir/Ah) 1 ^ 2 ; J^ 2 d(f>sm 3 (f>(l + cos 2 <fi)~ 2 = 1/2 
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